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Abstract. Fortran subroutines to calculate helicity amplitudes with massive spin-2 particles (massive gravi- 
tons), which couple to the standard model particles via the energy momentum tensor, are added to the 
HELAS (HELicity Amplitude Subroutines) library. They are coded in such a way that arbitrary scattering 
amplitudes with one graviton production and its decays can be generated automatically by MadGraph 
and MadEvent, after slight modifications. All the codes have been tested carefully by making use of the 
invariance of the helicity amplitudes under the gauge and general coordinate transformations. 



1 Introduction 

The idea of extra space dimensions has attracted much 
attention in recent years, since it can give us a novel so- 
lution to the hierarchy problem, or an alternative expla- 
nation of the hierarchical difference between the Planck 
scale (Mpi ^ 10^^ GeV) and the electroweak scale {mw ^ 
10^ GeV). 

So far, there have been various extra dimension mod- 
els, which can be divided into two major classes accord- 
ing to the geometry of the background space-time mani- 
fold. The first one includes the ADD (Arkani-Hamed, Di- 
mopoulos, and Dvali) model Ij and its variants, which ex- 
tend the dimension of the total space-time to D = A + 6, 
with a factorizable metric and large size of the compact 
extra dimensions (^ 1/Mpi). The second one includes the 
5-dimensional RS (Randall and Sundrum) model ^. and 
its variants, in which a warped metric is introduced along 
the 5-th dimension and the size of the extra dimension 
needs not to be much larger than the Planck length. 

In both classes of extra dimension models, there ap- 
pear Kaluza-Klein (KK) towers of massive spin-2 gravi- 
tons, which can interact with the standard model (SM) 
fields. The effective interaction Lagrangian is given by [31 

a 
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where T'"^^'' is the n-th graviton KK modes, and A is the 
relevant coupling scale. In the ADD model we have 



(2) 



A = Mpi = Mpi/V87r - 2.4 x W GeV, 



where My\ is the 4-dimensional reduced Planck scale, and 
in the RS model 



A = g-'^'^'^'^Mpi 



(3) 



is at the electroweak scale, where A; is a scale of order of 
the Planck scale and Tq is the compactification radius. 

In Eq. ([T]), T^i, is the energy- momentum tensor of the 
SM fields, 



= X\ 



H 



X. 



w 



(4) 



where 5^'' is the metric and r]^^ = diag(l, —1, —1, —1) is 
the Minkowski value, and each energy-momentum tensor 
is: 
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with e — gyy sin 9^ ~ 9z ^i^i '^^^ projection 
operator P^, = 5(1 ^ 75)- Here, the covariant derivative is 



Dp = d,, + igsT''A;^ + ieQfA,, 



(6) 



Note that the derivative couphngs of the W bosons are 
written exphcitly in Eq. (j5bp . Each field-strength tensor 
for the gauge bosons is 

f;, = dp a: - d^Al - gsF'^AlAl, 

Fp, = dp A, - d,Ap + ie[W+W- - in ^ i^)] , 

Zpu = dpZ^ - d„Zp + ig^ cosO^ \WpW^ - (a* ^ v)] , 



W^A, + cos e^w^z, -i^i^ u)] . 

(7) 

Notice that as in the standard HELAS package [S], we use 
the unitary gauge for the massive vector-boson propaga- 
tors and the Feynman gauge for the massless ones. 

In this paper, we present new HELAS subroutines [5] for 
the massive gravitons and their interactions based on the 
effective Lagrangian of Eq. ([T]), and implement them into 
MadGraph/MadEvent (MG/ME) [6 l [7 l [8]FI 



^ The Fortran code for simulations of the massive gravitons 
is available on the web [9]. 
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(5c) Table 1. List of the new vertex subroutines in HELAS system. 



The paper is organized as follows: In Sect. [2] we give 
the new HELAS subroutines. In Sect. [3] we describe how to 
implement amplitudes with a massive spin-2 graviton into 
MG/ME. In Sect, m we give sample numerical results. Sect. [5] 
contains a brief conclusion. 



2 HELAS subroutines for spin-2 particles 

In this section, we list the contents of all the new HELAS 
subroutines that are needed to evaluate massive spin-2 
gravtion production at hadron colliders in association with 
quark and gluon jets, and its decays into a pair of all 
the SM particles, or into arbitrary numbers of quarks and 
gluons. 

To begin with, in Sect. [2Tl1 the subroutine to compute 
the external lines for a spin-2 tensor particle is presented. 
Next, in the following subsections. 12.21 to 12. 7i the subrou- 
tines to compute the interactions of the tensor boson with 
the SM particles are explained. The new vertex subrou- 
tines are listed in Table [TJ Those subroutines which we do 
not present in this paper are FFST, SSST and SSSST type 
vertices for the interactions with scalars; VVVT and VVVVT 
for the interactions with the electroweak gauge bosons. 
These contribute e.g. to the graviton decays into three or 
more weakly interacting particles. We present the effective 
Lagrangian of Eqs. (5a) to (5f) for completeness sake. Fi- 
nally, we briefly mention how we test our new subroutines 
in Sect.ET 



2.1 Tensor wavefunction 

2.1.1 TXXXXX 

This subroutine computes the spin-2 Tensor particle wave- 
function; namely, e'^'^{p,X) and e^''(p. A)*, in terms of its 
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four-momentum p and helicity A, and should be called as 

CALL TXXXXX(P, TMASS, NHEL, NST , TC) 

The input P(0:3) is a real four-dimensional array which 
contains the four- momentum of the tensor boson, TMASS 
is its mass, NHEL {— ±2, ±1, 0) specifies its helicity A, and 
NST specifies whether the boson is in the final state (NST = 
1) or in the initial state (NST = -1). The output TC(18) 
is a complex f 8-dimensional array, among which the first 
f 6 components contain the wavefunction as 



TC(4/i 



namely 



where 



1^ + 1) 


= T(/i + 


TC( 1) 


-T(l,l), 


TC( 2) 


= T(1,2), 


TC( 3) 


= T(1,3), 


TC( 4) 


= T(1,4), 


TC(16) 


= T(4,4) 











(8) 



(9) 



and the last two contain the flowing-out four-momentum 

(TC(17), TC(18)) = NST (P(0) + iP(3), P(l) + iP(2)). 

(10) 

The helicity states of the tensor boson can be expressed 



with 



2.2 SST vertex 

The SST vertices are obtained from the interaction La- 
grangian among the tensor and two scalar bosons: 

£ssT = GTT^"^* [ - T^^,{{dPS*y{dpS*) - miss*] 

-f {dpS*y{d,s*) + [d^sy [d^s*) 

(16a) 

for the complex scalar field, or 

£ssT = GIT^''* [ - fl^.u{\^''S*^pS* - ^5*5*} 

+ d^,S*d^S*'\ (16b) 

for the real scalar field. Here, 

GT = GTS = -l/yl (17) 
is the coupling constant. 



e^''(p,±2) = e^(p,±)e''(p,±), 

e^''(p, ±1) = -i= [e^P, ±) e^b, 0) + e''(p, 0) e^ij,, ±)] , 
e^^ip, 0) = ^ [e^P, +) e'^b, -) + e^P, -) e'^b, +) 



2e''(p,0)e''(p,0)], 



(11) 



by using the vector boson wavefunctions (p. A) that obey 
the relation 



J_e''(p,A)=e^(p,A-l), 



(12) 



where J_ = Jx^iJy is the Jz lowering operator. The spin- 
1 vector wavefunction in the HELAS convention [5] satisfies 
this relation, and hence we simply use the HELAS code to 
obtain the tensor wavefunction. These tensor wavefunc- 
tions are traceless, transverse, orthogonal, and symmetric, 

e{p, = 0, p^e^^p, A) = p.e'"'(p. A) = 0, 
e^'^(p,A)e^.(p,A')*=<5Av, e^^- (p. A) = e'''' (p. A) , (13) 

and the completeness relation is 

+2 



2.2.1 SSTXXX 

This subroutine computes the amplitude of the SST vertex 
from two Scalar boson wavefunctions and a Tensor boson 
wavefunction, and should be called as 

CALL SSTXXX(S1, S2, TC, GT, SMASS , VERTEX) 

The inputs SI (3) and S2 (3) are complex three-dimensional 
arrays which contain the wavefunctions of the Scalar bosons, 
Sl(l) and S2(l), and their four- momenta as 

= (5ReSl(2),5ReSl(3),3TOSl(3),5mSl(2)), 
pt^ = (5ReS2(2),5ReS2(3),3TOS2(3),5mS2(2)). 

The input TC (18) is a complex 18-dimensional array which 
contains the wavefunction of the Tensor boson, and its 
four- momentum; see the TXXXXX subroutine in Sect. I2.1?l] 
GT is the coupling constant in Eq. (fTT)) in units of GeV""'^, 
and SMASS is the scalar boson mass ms in GcV units. Note 
that all the coupling constants in the latest HELAS library 
are defined as a complex number except for the 3-point 
and 4-point vector boson couplings. The output VERTEX is 
a complex number in units of GeV: 



J2 e''"(p,A)e"^(p,A)* - B^'^'^^ip) (14) VERTEX = GTT^''(m|?7^, - C^,,p,p^p^) Sl(l) S2(l), 



A=-2 



(18) 
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where we used the notation 

T'^'^ = T(/i+l,j/+l)=TC(4M + i^+l) (19) 
in Eqs. ^ and ®, and 

(20) 

2.2.2 HSTXXX 

This subroutine computes an off-shell scalar current H 
made from the interactions of a Scalar boson and a Tensor 
boson by the SST vertex, and should be called as 

CALL HSTXXX(TC, SC, GT, SMASS, SWIDTH , HST) 

The inputs TC(18) and SC(3) are the wavcfmictions and 
momenta of the Tensor and Scalar bosons, respectively, 
and SWIDTH is the scalar boson width Fs- The output 
HST (3) gives the off-shell scalar current multiplied by the 
scalar boson propagator, which is expressed as a complex 
three-dimensional array: 



HST(l) = iGT 



rpfll/ 



(f- - mg -I- imsTs 

X (m|77^, + C^^.p^pPq^) SC(1), (21) 



and 



HST(2) ^ TC(17) + SC(2), 
HST(3) = TC(18) + SC(3). 



(22) 
(23) 



Here the momenta p and q are 



= (3?eSC(2),5ReSC(3),3TOSC(3),3mSC(2)), 

= (3fJeHST(2),5ReHST(3),3mHST(3),3mHST(2)). 

2.2.3 USSXXX 

This subroutine computes an off-shell tensor current U 
made from two flowing-out Scalar bosons by the SST ver- 
tex, and should be called as 

CALL USSXXX(S1, S2, GT, SMASS, TMASS, TWIDTH , USS) 

The inputs TMASS and TWIDTH are the tensor boson mass 
and width, ttit and Ft- The output USS (18) gives the off- 
shell tensor current multiplied by the tensor boson prop- 
agator, which is expressed as a complex 18-dimensional 
array: 



— m'^p + irriTFT 
X {mlri,,, - C^.,p.KpD S1(1) S2(1), (24) 



where we used the notation 

T"'3 = T(a-fl,/3+l), 



whose components are assigned into the first 16 compo- 
nent of USS as in Eq. ([5]), and 



USS(17) = Sl(2) + S2(2), 
USS(18) Sl(3) + S2(3). 



(26) 
(27) 



Here, pi and P2 are the momenta of the outgoing scalars, 
and q is that of the off-shell tensor boson given in Eqs. (pS)) 
and (1271) as 



(25) 



q^' = (3fieUSS(17),3?eUSS(18),3mUSS(18),3mUSS(17)). 

Although the effective Lagrangian of Eq. ^ does not 
dictates the off-shell behavior of the gravitons, we allow 
gravitons to propagate just once in the total amplitude 
where there are no external gravitons in the initial or fi- 
nal states. This is convenient when studying the correlated 
decays of the graviton production and its subsequent de- 
cays. 

We may also note that the order of the GT couplings 
should be restricted to 1 when there is an external gravi- 
ton, and 2 when a graviton is exchanged among the SM 
particles. 

Before turning to the FFT vertex, it should be noticed 
here that the conventional factors of i in the vertices and 
those in the propagators are both included in the off-shell 
wavefunctions, such as Eqs. (f2T|) and (f24l) above, according 
to the HELAS convention. The HELAS amplitude, obtained 
by the vertices, such as Eq. (fT8| . gives the contribution 
to the T matrix clement without the factor of i. See more 
details in the HELAS manual [5j. 



2.3 FFT vertex 

The FFT vertices are obtained from the interaction La- 
grangian among the tensor boson and two fcrmions: 

+ {\hl,d.I - \d^{hluf) + (m v)}] , (28) 
where the coupling constant is 

GT = GTF = -l/(4yl). (29) 



2.3.1 lOTXXX 

This subroutine computes the amplitude of the FFT vertex 
from a flowing-In fermion spinor, a flowing-Out fermion 
spinor and a Tensor boson wavefunction, and should be 
called as 

CALL IOTXXX(FI,FO,TC, GT, FMASS , VERTEX) 

The inputs FI(6) and F0(6) are complex six-dimensional 
arrays which contain the wavefunctions of the flowing-In 
and flowing-Out Fermions, and their four-momenta as 

= (3?eFl(5),g?eFl(6),3TOFl(6),3mFl(5)), 
pt^ = (3fJeFD(5),3fJeF0(6),3TOFD(6),3mF0(5)). 
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The input GT is the coupUng constant in Eq. ((29|) . and 
FMASS is the fermion mass mp. What we compute here is 

VERTEX = GT (FD) [T^^" { - ri^,„{^^+ ^3 - "^mp) 

+ 7Mbi+P2). + (M^'^)}](Fl), (30) 

where we use the notation 



(FI) = 



'FI(1)~ 

FI(2) 

FI(3) J ' 
.FI(4), 



(31) 



(FO) = (FD(1), F0(2), FD(3), FD(4)). (32) 



2.3.2 FTIXXX 

This subroutine computes an ofF-shell Fermion wavefunc- 
tion made from the interactions of a Tensor boson and a 
flowing-In fermion by the FFT vertex, and should be called 
as 

CALL FTIXXX(FI,TC,GT, FMASS, FWIDTH , FTI) 

The output FTI (6) gives the off-shell fermion wavefunc- 
tion multiplied by the fermion propagator, which is ex- 
pressed as a complex six-dimensional array: 

(FTI) - ,rT '-(^+"^^) T'^- 

+ -f^ip + q), + ifi^iy)]{Fl), (33) 
where we use the notation 



and 



(FTI(1)~ 
FTI (4), 



FTI(5) FI(5) - TC(17), 
FTI(6) = FI(6) - TC(18). 



(34) 



(35) 
(36) 



Here the momenta p and q are 



= (3fJeFl(5),5ReFI(6),3TOFl(6),5mFl(5)), 
qf" = (3fieFTl(5),3?eFTl(6),3mFTl(6),3mFTl(5)). 

2.3.3 FTDXXX 

This subroutine computes an off-shell Fermion wavefunc- 
tion made from the interactions of a Tensor boson and 
a flowing-Dut fermion by the FFT vertex, and should be 
called as 

CALL FTOXXX(FO, TC, GT, FMASS, FWIDTH , FTD) 



The output FTQ(6) is a complex six-dimensional array: 

(FTD) = iGT (FD) T^'' 

X [- 4~ 2mF) +lt_i[p + q)v + (m ^ ^)\ 

i{i + mF) 



q^ — rn?p + impFp ' 



(37) 



where we use the notation 

(FTO) = (FTO(l), FT0(2), FT0(3), FT0(4)), (38) 

and 



FT0(5) = F0(5) + TC(17), 
FT0(6) = F0(6) + TC(18). 



(39) 
(40) 



Here the momenta p and q are 

p'' = (5ReF0(5),3?eF0(6),3mF0(6),3TOFD(5)), 
= (3fieFT0(5),3f?eFT0(6),3mFT0(6),3mFT0(5)). 



2.3.4 UIOXXX 

This subroutine computes the bi-spinor tensor current U 
made from flowing-In and flowing-Dut fermions by the FFT 
vertex, and should be called as 

CALL UIDXXX(FI, FQ, GT, FMASS, TMASS, TWIDTH , UID) 

The output UI0(18) is a complex 18-dimensional array: 



q^ — mi, -f imTTT 



X (FQ)[-r/^^(^i+^2-2mF) 

+ 7m(Pi+P2). + (m^'^)](FI) (41) 

for the first 16 components of UID, and 



UI0(17) = -FI(5) +FD(5), 
UI0(18) = -FI(6) +FD(6). 



(42) 
(43) 



Here, pi and p2 are the momenta of the flowing-in and 
flowing-out fermions, respectively, and q is that of the ten- 
sor particle. 



2.4 VVT vertex 

The VVT vertices are obtained from the interaction La- 
grangian among the tensor and two vector bosons: 



■■} 



(44a) 



6 



K. Hagiwara, J. Kanzaki, Q. Li and K. Mawatari: HELAS and MadGraph/MadEvent with spin-2 particles 



for the complex vector bosons, or W bosons, and 



2.4.2 JVTXXX 



GT '^^^ 



^|i(aPT/" - d''vp*){dpv: ~ d,v;) 

[d^yp* ~ dpv;){d,v; - d^v:) + mlv;v: 
C'Ti,,(d''d''v:v; + ]^d''v;d''v: 



(44b) 



for the real ones, or gluons, photons, and Z bosons. Here 
the coupling constant is 



GT = GTV : 



A/A. 



(45) 



The ^ terms are the gauge-fixing terms, which vanish for 
massive vector bosons in the unitary gauge. For massless 
vector bosons we take ^ = 1 in the Feynman gauge. 



2.4.1 VVTXXX 

This subroutine computes the amplitude of the VVT vertex 
from two Vector boson polarization vectors and a Tensor 
wavefunction, and should be called as 

CALL VVTXXX(V1,V2,TC, GT,VMASS , VERTEX) 

The inputs VI (6) and V2(6) are complex six-dimensional 
arrays which contain the Vector boson wavefunctions, and 
their momenta as 

= (5ReVl(5),3fieVl(6),3mVl(6),3TOVl(5)), 
= (5ReV2(5),5ReV2(6),3mV2(6),3TOV2(5)). 

The input GT is the couphng constant in Eq. and 
VMASS is the vector boson mass mv- What we compute 
here is 

VERTEX = GTT'"'[(m^ + Pi ■ P2) C^^^p„ + I^^^.p^ (^1,^2) 
+ i-^E^,,,p„{Vi.V2)\VlV^. (46) 



where we use the notation 



Fi^ = Vl(Ai+l), 
K,^ = V2(a*+1), 



(47) 
(48) 



and 



Df_,„,pa{pi,P2) = r]f^iyPl<yP2p - [rjp.crPlvP2p + r]iLpPlaP2u 

- 'np<yPipP2u + (a^ ^ v)] , (49) 

Ef^u,pa{Pl,P2) = Vp-uiPlpPla + P2pP2a + PlpP2a) 

- [Vi^aPlfiPlp + VupP2pP2a + (A* ^ ^)\ ■ 

(50) 



This subroutine computes an off-shell vector current J 
made from the interactions of a Vector boson and a Tensor 
boson by the VVT vertex, and should be called as 

CALL JVTXXX(VC, TC, GT, VMASS, VWIDTH , JVT) 

The input VC(6) is the wavefunction and momentum of 
the Vector boson. The output JVT(6) gives the ofF-shell 
vector current multiplied by the vector boson propagator, 
which is expressed as a complex six-dimensional array: 



JVT(a + 1) = iGT 



q- — rriy + irnvry 



X T^'" [{ml -p-q) Cp,,,pa + D^,,,p„{p, -g)] V^' 

(51) 

for the massive vector boson, or 

JVT(a +1)^ iQT ^ ry^^T^" [- ip ■ q) Cp,,pa 

+ D^,,p, {p, -q) + E^,^p„ (p, -q)] VP (52) 
for the massless vector boson, and 



JVT(5) = VC(5) -HTC(17), 
JVT(6) = VC(6) -KTC(18). 



(53) 
(54) 



Here the momenta p and q are 

= (5ReVC(5),3fieVC(6),3mVC(6),3TOVC(5)), 
q^" = (3fieJVT(5),5ReJVT(6),3TOJVT(6),3mJVT(5)). 



2.4.3 UVVXXX 

This subroutine computes an off-shell tensor current U 
made from two flowing-out Vector bosons by the VVT ver- 
tex, and should be called as 

CALL UVVXXX(V1, V2, GT, VMASS, TMASS, TWIDTH , UVV) 

The output UVV (18) is a complex 18-dimensional array: 

■ f>pu,af3 
na(3 ■ — '-O 



yap ^ -Q^ 



q^ — m3p + imTFT 
X [{m?y +pi ■ P2) Cp^.pa + Dp^,p^(pi,p2) 
+ r'E^,^p^{p,,P2)]VfV2^ (55) 



for the first 16 components of UVV, and 

UVV(17) = Vl(5) + V2(5), 
UVV(18) = Vl(6) + V2(6). 



(56) 
(57) 



Here, pi and p2 are the momenta of the outgoing vector 
bosons, and q is that of the tensor boson. 
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2.5 FFVT vertex 



2.5.3 FVTOXX 



The FFVT vertices are obtained from the interaction La- 
grangian among the tensor boson, vector boson and two 
fermions: 

£ffvt = -2GTT^'^*(?7^^.77p^ - rj^^rj^p) 

X /7- [GC(1)Pl + GC(2)Pfl] / VP* (58) 

with the chiral-projection operator Pr l — 5(1 ±75)- The 
couphng constant GT is 



GT = GTFV = -1/(2/1), 



(59) 



and GC(1) and GC(2) are the relevant FFV left and right 
couphng constants. The hst of the coupling constants is 
shown in Table [21 For instance, in the case of the interac- 
tion with a gluon, the FFV couplings are 



GC(1) = GC(2) = -g. 



(60) 



where the sign of the coupling is fixed by the HELAS con- 
vention [5j. 



This subroutine computes an off-shell fermion wavefunc- 
tion by the FFVT vertex, and should be called as 

CALL FVTOXX(FO, VC, TC, GC, GT, FMASS, FWIDTH , FVTO) 

What we compute here is 

(FVTO) = -^GT {FO)TP''{tj^,tJp^ - Cp,,p„) 

X [GC(1)P^ + GC(2)P;,] K^+rnp) 

— nip + irapl p 
(65) 



for the first 4 components of FVTO (6), and 

FVT0(5) = F0(5) + VC(5) + TC(17) 
FVT0(6) = F0(6) + VC(6) + TC(18) 



(66) 
(67) 



Here the momentum q is 

= (5ReFVT0(5),5ReFVT0(6),3mFVT0(6),9mFVT0(5)). 



2.5.1 IDVTXX 

This subroutine computes the amplitude of the FFVT ver- 
tex from a flowing-In fermion spinor, a flowing-Out fermion 
spinor, a Vector boson polarization vector and a Tensor 
boson wavefunction, and should be called as 

CALL IOVTXX(FI, FO, VC, TC, GC, GT , VERTEX) 

What we compute here is 

VERTEX - GT (FO) [r'^''(?7^^77p^ - Cp^,p^) 

X 1/V(GC(1)Pl + GC(2)Pj^)] (FX). (61) 

2.5.2 FVTIXX 

This subroutine computes an off-shell fermion wavefunc- 
tion by the FFVT vertex, and should be called as 

CALL FVTIXX(FI, VC, TC, GC, GT, FMASS, FWIDTH , FVTl) 

What we compute here is 

(FVTI) = -^GT ^ ^(^+"^^) T^-(ryp.77p.-Cp.,p.) 
q'' — nip + iniF-L F 

X ypj^ [GC(1)Pl + GC(2)Pi^] (Fl) (62) 

for the first 4 components of FVTI (6), and 

FVTI(5) FI(5) - VC(5) - TC(17), (63) 
FVTI(6) = FI(6) - VC(6) - TC(18). (64) 

Here the momentum q is 

= (5ReFVTl(5),3fieFVTl(6),3mFVTl(6),3mFVTl(5)). 



2.5.4 JIOTXX 

This subroutine computes an off-shell vector current by 
the FFVT vertex, and should be called as 

CALL JIOTXX(FI, FO, TC, GC, GT, VMASS, VWIDTH , JIOT) 

What we compute here is 



JIOT(q;-F 1) = -iGT- 



mi 



q^ — my + irnvTy 

X (FO) [7'^{GC(1)Pl + GC(2)Pfl}] (Fl) (68) 
for the massive vector boson, or 

JIOT(a + 1) = -iGT TjP" Tf^'iTjp.rip, - Cp,,pa) 

X (FO) [7"{GC(1)Pl + GC(2)P,^}] (Fl) (69) 

for the massless vector boson, and 

JI0T(5) = -FI(5) + F0(5) + TC(17), (70) 
JI0T(6) = -FI(6) + F0(6) + TC(18). (71) 

Here q is the momentum of the vector boson. 
2.5.5 UIOVXX 

This subroutine computes an off-shell tensor current by 
the FFVT vertex, and should be called as 

CALL UIOVXX(FI, FO, VC, GC, GT, TMASS, TWIDTH , UIOV) 
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What we compute here is 



— + irriTrT 
X (FO) [1/''7"{GC(1)Pl + GC(2)Pi?,}] (Fl) (72) 
for the first 16 components of UIDV(18), and 

UI0V(17) = -FI(5) + F0(5) + VC(5), (73) 
UIDV(18) = -FI(6) + F0(6) + VC(6). (74) 

Here q is the momentum of the tensor boson. 
2.6 VVVT vertex 

The VVVT vertices are obtained from the interaction La- 
grangian among the tensor and three vector bosons: 



1 



rb,p-irc,a 



(75) 

with the structure constant /"'"^ and the couphng con- 
stant, as in Eq. (|45)) . 



GT = GTV : 



-l/A. 



(76) 



In this paper we concentrate on the interactions with 
gluons for the VVVT vertex, so in this case GC is the strong 
couphng constant, 



GC 



9s, 



(77) 



and /"^'^ is the structure constants of the group SU{i), 
which can be handled by the MG automatically. As in the 
original subroutines for the VVV vertex [5], the following 
subroutines, VVVTXX, JVVTXX, and UVVVXX, can be used for 
the electroweak gauge bosons without any modifications. 



2.6.2 JVVTXX 

This subroutine computes an off-shell vector current by 
the VVVT vertex, and should be called as 

CALL JVVTXX(VA, VB, TC, GC, GT, VMASS, VWIDTH , JVVT) 

What we compute here is 

JVVT (a + 1) 



-iGT GC 



g2 - + imyrv 



V 



{pbp + qp) 

(80) 

for the massive vector boson, or 

JVVT(a+ 1) - -iGTGC-j77^"T'"'[C^,,p,(paA - Pbx) 

+ C^i/.crA (Pbp + Qp) + C^y^\p{-qa - Paa) 

+ F^,,p,x{Pa,Pb,~q)\VPV,% (81) 

for the massless vector boson, and 

JVVT(5) = VA(5) + VB(5) + TC(17), (82) 
JVVT(6) = VA(6) -h VB(6) -|-TC(18). (83) 

Here the momenta Pa , Pb and q are 

= (3fieVA(5),3fieVA(6),QTOVA(6),3mVA(5)), 
= (3fieVB(5),3?eVB(6),9mVB(6),3mVB(5)), 
= (5ReJVVT(5),5ReJVVT(6),3TOJVVT(6),3TOJVVT(5)). 

Note that the off-shell gluon JVVT has the color c. 



2.6.1 VVVTXX 

This subroutine computes the amplitude of the VVVT ver- 
tex from three Vector boson polarization vectors and a 
Tensor boson wavefunction, and should be called as 

CALL VVVTXXX(VA, VB, VC, TC, GC, GT , VERTEX) 

What we compute here is 

VERTEX = - GT GC T'"' [C^,,p, [pax - Pbx) 

+ C pv,ax(pbp Pep) + C p.v,Xp(pca Paa) 

+ F^,,p,x{pa,Pb,Pc)]VPV,^V,^ (78) 

with 

Ff,u,p<7x{Pa,Pb,Pc) = VppV<7x{Pb-Pc)iy+VpaVpx{Pc-Pa)u 
+ Vp-XVpcriPa ~Pb)u + (Ai ^ 1^)- (79) 

Here, the vector bosons (gluons in this paper) VA, VB, and 
VC have the momentum pa, pt, and pc, and the color a, b, 
and c, respectively. 



2.6.3 UVVVXX 

This subroutine computes an off-shell tensor current by 
the VVVT vertex, and should be called as 

CALL UVVVXX(VA, VB, VC, GC, GT, TMASS, TWIDTH , UVVV) 

What we compute here is 

T"^ = -iGTGC '^"'^'"'^ ^ [C^.,pa(PaA - Pbx) 

+ C^i,^„x{Pbp - Pep) + Cp^^^xpip ca Paa) 
+ Fp,,p^x{pa,Pb,Pc)] V^V.^V,^ (84) 

for the first 16 components of UVVV (18), and 

UVVV(17) = VA(5) + VB(5) +VC(5), (85) 
UVVV(18) = VA(6) + VB(6) +VC(6). (86) 

Here Pa , Pb and Pc are the momenta of the outgoing vector 
bosons, and q is that of the tensor boson. 
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2.7 VVVVT vertex 



and 



The VVVVT vertices are obtained from the interaction La- 
grangian among the tensor and four vector bosons: 



VVVVT — — bi IjU / / 
1 



i2 4^abe ^cde rpfj,i/^ 

r]^^v''^p*v'''''*v;;*vt - v'''P*v;^*v^"v;^* 



(87) 



with the couphng constants GT GTV in Eq. (gS]) and GC 
= (7s for the interactions with gluons. 

We should note that the 5-point vertex cannot be gen- 
erated by MG, and thus we must add the following sub- 
routines, GGGGTX, JGGGTX, or UGGGGX, to the amphtudes 
which have the corresponding color structures by hand. 



2.7.1 GGGGTX 

This subroutine computes the portion of the amplitude of 
the VVVVT vertex from four Gluon polarization vectors and 
a Tensor boson wavefunction corresponding to the color 
structure f°-be jcde^ should be called as 

CALL GGGGTX(VA, VB, VC, VD, TC, GC, GT , VERTEX) 

The output is 

VERTEX = -GT GC^ T^'' G^,,pXaS VPV{^V^V! (88) 

with 

+ [Vtip'nvS'riXa- + TlpWvallpS ~ VppVvcrllXS 

- VfiXVvsrip'y + (m ^ ^)\ ■ (89) 

To obtain the complete amplitude, this subroutine must 
be called three times (once for each color structure) with 
the following permutations: 

CALL GGGGTX(VA, VB, VC, VD, TC, GC, GT , VERTEXl) 
CALL GGGGTX(VA, VC, VD, VB, TC, GC, GT , VERTEX2) 
CALL GGGGTX(VA, VD, VB, VC, TC, GC, GT , VERTEX3) 

corresponding to the color structure f^bej^cde^ jace jdbe ^ 
and f^dejbce^ respectively. 

2.7.2 JGGGTX 

This subroutine computes the portion of the off-shell gluon 
current by the VVVVT vertex, corresponding to the color 
structure f^bejcde^ should be called as 

CALL JGGGTX(VA, VB, VC, TC, GC, GT , JGGGT) 

What we compute here is 

JGGGT(a + 1) = -zGTGC^ ^ 77*" T^'' G,,,,pXa5 V^V^^V^, 

(90) 



JGGGT(5) = VA(5) + VB(5) + VC(5) + TC(17), (91) 
JGGGT(6) = VA(6) + VB(6) + VC(6) + TC(18). (92) 

Note that the off-shell gluon JGGGT has the color d and 
the momentum q. 

2.7.3 UGGGGX 

This subroutine computes the portion of the off-shell ten- 
sor current by the VVVVT vertex, corresponding to the color 
structure jo^be jcde ^ should be called as 

CALL UGGGGX(VA, VB, VC, VD, GC, GT, TMASS, TWIDTH , UGGGG) 

What we compute here is 



T"^ = -iGTGC^- 



— + irriTrT 



(93) 



for the first 16 components of UGGGG (18), and 



UGGGG(17) = VA(5) + VB(5) + VC(5) + VD(5), (94) 
UGGGG(18) = VA(6) + VB(6) + VC(6) + VD(6). (95) 

Here q is the momentum of the tensor boson. 



2.8 Checking for the new HELAS subroutines 

The new HELAS subroutines are tested by using the QCD 
gauge invariance and the general coordinate transforma- 
tion invariance of the helicity amplitudes. In particular, 
we use the following processes; 

qq gT for lOTXXX, FTIXXX, FTDXXX, lOVTXX, 

gg gT for VVTXXX, JVTXXX, VVVTXX, 

QQ qq{99)T for FVTIXX, FVTDXX, JIDTXX, 
gg ggT for JVVTXX, GGGGTX, 

QQigg) ^ gg for UIOXXX, UVVXXX. 

More explicitly, we express the helicity amplitudes of the 
above processes as 



(96) 



with an external tensor and a gluon wavefunction. The 
identities, 

Pt/PtT^^up e''{pg,Xg)* = Pt/PtT,iup eP{pg,Xg)* = 

(97) 

for the general coordinate transformation symmetry and 

P^/P°7^M^P = (98) 

for the SU(3) gauge invariance, test all the above subrou- 
tines thoroughly. 
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The subroutines UIDVXX, UVVVXX, JGGGTX and UGGGGX 
have been checked in such a way that the above subrou- 
tines are used rewriting the hehcity amphtudes of the pro- 
cesses 

qq^T ^ qqgiggg), 
gg ggT ggZZ, 

99 99g9T- 

We test the agreement of the amphtudes between the ex- 
pressions that uses the above subroutines and those with- 
out them for all helicity combinations and at arbitrary 
Lorentz frame. 



3 MadGraph/MadEvent implementing for 
spin-2 gravitons 

In this section, we would like to describe how we imple- 
ment spin-2 gravitons into MG/ME. 

First, using the User Model framework in MG 8 , we 
make our new model directories for both the ADD and RS 
models, including the massive gravitons (particles . dat) 
and their interactions with the SM particles (couplings . f 
and interactions .dat); see also TableO Then we insert 
all the new HELAS subroutines for spin-2 tensor bosons into 
the HELAS library in MG. Since the present MG does not han- 
dle external spin-2 particles, we further modify the codes 
in MG to tell it how to generate the SST, FFT and FFVT type 
of vertices and helicity amplitudes (for VVT and VVVT type, 
it has already been done for the Higgs effective field theory 
(HEFT) model), and how to deal with the helicity of the 
spin-2 tensor bosons when they are external. Moreover, 
since MG can only generate Feynman diagrams with up 
to 4-point vertices, the amplitudes and their HELAS codes 
with the 5-point vertex, GGGGTX, JGGGTX, or UGGGGX, have 
been added by hand; see more details in Sect. 12.7] 

Finally, we note that since in the ADD model the gravi- 
tons are densely populated and we should sum over their 
contributions by modifying the phase space integration 
in ME. In the ADD model, the spectrum of KK graviton 
modes can be treated as continuous for (5 < 6 [I], and the 
mass density function is given by [3] 

27r'^/2 

p{m) ^ Ss-ij^m^-'^ with Sg-i = p(^s/2) ' ^^^^ 

where Mg is the ADD model effective scale. Thus, we mod- 
ify the phase space generating codes in ME to add one more 
random number for graviton mass generating and imple- 
ment the above graviton mass integration. 

4 Sample results 

In this section, we present some sample numerical results 
for the graviton plus mono-jet and di-jet productions at 
the LHC, using the new HELAS subroutines and the mod- 
ified MG/ME. 
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Table 2. List of the coupling constants for each vertex. All 
the particles and the coupling constants are written in the MG 
notation, y stands for a massive graviton, f represents all pos- 
sible fermions, and v is the SM gauge bosons (g, a, z, w). GC 
is a SM coupling constant, while GT is a non-renormalizable 
coupling constant defined in each subroutine in Sect. [D 



We use the following jet definition criteria 



P^>20GeV, \vj\<5, Rjj = \/ Aq^ + Z\02 > o.6, 

(100) 

where 77 is the pseudorapidity of the jets and <^ is the 
azimuthal angle around the beam direction, and further 
require 

-^32 < 0, \V3,~Vj2 \ > 4.2. (101) 

The CTEQ6L1 parton distribution functions [11] are em- 
ployed with the factorization scale chosen as /i/ = min(Py ) 
of the jets which satisfy the above cuts. For the QCD 
coupling, we set it as the geometric mean value, = 

Fig. [1] shows the P^"^*' distributions for the graviton 
productions with the mono-jet (via the gg, qg and qq chan- 
nels) and with the di-jet (via the gg, qg and qQ channels) 
in the ADD model at the LHC, with A = 5 TeV and 
(5 = 4. Note that the unitarity criterion Mt„ < /I is used. 
See more details in Ref. [TO] . 

Fig. [2] shows the (pjj distributions for the first KK 
graviton productions with two jets in the RS model at 
the LHC, via the gg, qg and qQ channels. Here we set 
yl = 4 TeV and Mt^ = 1 TeV. The total cross sections 
are 0.74, 3.74 and 4.33 pb for the qQ, gg and qg channels, 
respectively. 
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Fig. 1. P™"''^ dependence of the total cross sections for the 
graviton productions with 1-jet (via the gg, qg and qq channels) 
and 2-jets (via the gg, qg and qQ channels) in the ADD model 
at the LHC, where q,Q = u, d, s, c. 
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Fig. 2. Distributions of the azimuthal angle separation be- 
tween the two jets for the first KK graviton plus di-jet pro- 
ductions in the RS model at the LHC, via the gg, qg and qQ 
channels (g, Q — u,d, s, c). 



5 Summary 

In this paper, we have added new HELAS subroutines to 
calculate helicity amplitudes with massive spin-2 particles 
(massive gravitons) to the HELAS library. They are coded 
in such a way that arbitrary scattering amplitudes with 
one graviton production and its decays can be generated 
automatically by MG and ME, after slight modifications. All 
the codes have been tested carefully by making use of the 
invariance of the helicity amplitudes under the gauge and 
general coordinate transformations. 



